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Abstract 

In this note, we present a detection criterion for exo-planets to be used with the space mission 
COROT. This criterion is based on the transit method that suggests to look for stars dimming 
caused by partial occultations by planetary companions. When at least three transits are observed, 
we show that a cross-correlation technique can yield a detection threshold, thus enabling to evaluate 
the number of possible detections assuming a model for the stellar population in the Galaxy. 
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1 Overview of COROT 

1.1 Introduction 

Exo-planets detection has now become a very active field in astrophysics. In august 2000, about fifty 
Jupiter-Hke planets have been discovered around nearby stars. The next challenging step is to find much 
less massive planets with the hope to detect life on some of them afterwards. COROT, a CNES space 
mission to be launched in 2004 [1], is currently a funded project capable of detecting planets with a 
radius close to li?© thanks to the "transit method" [2]. In this section, we briefiy describe the transit 
method and the COROT satellite. In section 2, we present a detection criterion for multi-transits, and 
in section 3, we discuss the number of possible dectections assuming a model for the stellar population 
in the Galaxy. The specific case of mono-transits will be considered in a forthcoming article. 

1.2 The transit method 

The transit method for searching extrasolar planets is based on the idea that if a planet crosses the disk 
of its parent star, it would result in a dimming of the star's observed light [3-4]. The expected amplitude 
of the relative stellar fiux dimming is: 



where F is the flux, Rp the radius of the planet and i?* the radius of the parent star. We call impact 
parameter, expressed in stellar radii, the apparent height of the planet trajectory above the star's equator. 
For an impact parameter equal to 0.5, the transit duration is: 



tr = 



(2) 



where a is the orbital radius. This duration is respectively 11.2 and 25.7 hours for the Earth and Jupiter 
in the solar system, and 3 hours for 51 Peg b [5]. The geometrical probability that the orbital inclination 
on the sky is close enough to 90° to make a transit visible is: 



Pg equals to 0.5% for the Earth, 0.1% for Jupiter and 16% for 51 Peg b. This method has already been 
succesfully used from ground and space on the planetary companion of HD209458, previously detected 
by radial velocity techniques, e.g. [6]. 

1.3 Focusing on multi-transits 

COROT features a 27 cm telescope and four 2048 x 2048 CCDs, two of which are devoted to the exo- 
planets program. During its 2.5-year mission, the satellite will monitor 5 fields of 5000 to 12000 stars 
(11 < mv < 16.5, galactic latitude = 15 — 20°), each of them for 150 days, thus leading up to 60000 
lightcurves. The task will be then to look in the data for the signatures of planetary transits. In order 
to estimate the number of potential detections by COROT, we propose to consider here a rather basic 
method based on a cross-correlation technique. More accurate techniques are under study at the Labo- 
ratoire d' Astronomic Spatial in Marseille [7], or have been already suggested [8]. Here we will consider 
only multi-transits, which means that a given planet will have to transit at least 3 times in front of its 
parent star to be detected by this method. The sought after signal looks like a repeated dimming in the 
parent star lightcurve, at even time intervals corresponding to the planetary orbital period P. 

2 Principle of data processing 

2.1 Raw data averaging 

To increase the signal to noise ratio (SNR) and also to reduce the amount of data, we begin by averaging 
the raw data, i.e. samples of the total flux F every 16 minutes over 150 days, on the duration of a 
presumed transit tr (up to 15 hrs). This operation leads to a iV = 150 x 24/tr point vector in which 
a potential event is reduced to one point. It means that in the whole process of the transit search, all 
relevant values of tr should be tried. 

2.2 Transit signal modelling 

The transit signal s we are looking for can be modeled by the sum of a Dirac comb allk (k teeth, 
amplitude a) and a gaussian white noise b of standard deviation (76 : 
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{S is Kroenecker's symbol). The term b is the sum of several types of noises: 
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• the photon noise: aph = \/Nph {Np^ refers actually to the number of stellar photo-electrons detected 
by the detector during tr) 

• the electronical read-out noise: aro = 12 e~ .pixel~^ 

• the background noise (essentially zodiacal Hght): atg = 16 e~.pixel~^ (exposure time is 32 s) 

• the stellar irradiance variability noise depending on the averaging interval whose length is tr: 
o'st{ii') = Nph{tr) Ost{tr)ppm- For G and K stars, an estimated value of this quantity can be obtained 
through filtering of the SOHO-VIRGO data acquired on the Sun [9]: (Tst{i'r)ppm ~ 30 — 75 ppm for 
tr = 2— 15 hrs. 

• the noise introduced by random pointing error of the satellite is supposed to be perfectly corrected 
by onboard processing 

For simplicity, all those noises are supposed independant, white and gaussian noises. Consequently, we 
get: 

(^b = ^Nph + n{a^,, + alg) + a% (6) 
where n denotes the total number of read-out pixels during the transit. 



2.3 Detection with cross-correlation 



The idea is to compute cross-correlation products between the data and a Dirac comb (amplitude unity) 
which has the shape of a noise free multi-transit signal (see fig.l). For a given value of k, we get N/k 
products of this kind, denoted Ck- 
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Trying all values of k between 3 and 50 (3 < P < 50 days) means computing X^^^La N/ k = 'jN products 
(7 ~ 3.0). The presence of a transit must be somehow related to a high value of Ck, but the question is: 
where to draw the line? 



2.4 Statistical analysis 

The answer to this question lies in a statistical analysis of the problem. To begin with, let us consider the 
set of N/k products Ck for a given k. Ck can be seen as a random variable. Because we assumed that b 
is a gaussian noise with a null mean and a standard deviation a^, the probability law of Ck should also 
be gaussian with a standard deviation ac = Vkab/N. Its mean value should be equal to ak/N in case 
of a star actually showing transits or zero otherwise. Let pk be the probability of having Ck < /3/s crc in 
case of noise only (no transit): 

Pk = Pr{Ck < Pk(Tc} = ^ 

The probability that all the N/k values of Ck will remain inferior to (3kcrc is Ofe = Pk'"- This gives the 
level of confidence that the statistical noise would not generate a high value of Ck that would be mistaken 
for a transit (see fig. 2). If the statistics bear on all values of k {'yN cross-correlation products), the level 
of confidence becomes: 

50 50 , ^ 
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One can numerically show that (3 depends weakly on k provided a is close enough to one (no more than 
10% variations). That is why we will assume that (3 is common to all fc, in order to solve for /3, given the 
global level of confidence a: 

erf ^-^^ + 1 - 2oW = (10) 

For example, the transit of a planet orbiting a sun-like star at 0.05 AU (such as 51 Peg b) would last 
approximately 3 hrs if the line of sight belongs to the orbital plane. In this case N = 1200, and solving 
(|^ for a = 99.9 % leads io (3 = 5.0 (see fig.3). 

2.5 SNR detection criterion 

Assume that a global level of confidence a has been chosen and has yielded up a value of (3: a detection 
could be claimed with that confidence, if among the "fN cross-correlation products, one is greater than 
I3ac- As it is a signature of a transit, its value can be written as ak/N, so the detection criterion 
translates into the inequation: 

ak ^ VkPab 

Introducing the SNR on a single event (a single dimming of the Hghtcurve) defined by S/N — a /at,, we 
get: 

As expected, the required SNR increases with the level of confidence and decreases with the number of 
observed transits. To go on with the previous example, P — 4.1 days so A: = 36 and S/N > 0.7 ! As we 
see here, this cross-correlation technique is enough powerful to detect repetitive transits with lower than 
unity SNR on single events. 

2.6 Link with the planetary radius 

The amplitude of the relative dimming of the parent star Hghtcurve is connected to the radius of the 
transiting planet Rp by: 
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Thus (12) translates into: 



besides = = (14) 

F Nph NNph ^ ' 



which gives the minimum planetary radius that can be detected with a level of confidence a. In our 
example, if we assume for the star mv = 14 and given the characteristics of COROT, we compute: 
Rp > 1.5 i?©. 



3 Expected number of detections 
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3.1 Assumptions 

To have a realistic distribution of stars per spectral type and magnitude interval in the observed fields, we 
have used a model of the stellar population in the Galaxy developped at Besangon Observatory [10]. Only 
stars on the Main Sequence have been considered (type V stars). We will assume that every observed 
star has 100% probability of having a planet orbiting at the distance tested, and that orbits are circular. 
For every distance, we compute the minimum radius the planet should have for our criterion to be able 
to detect it with a false alarm level of 10~^. 

3.2 Algorithm 

Let a be the orbital of the planetary orbit. We define the reduced orbital radius by — a{Li,/ Lq) 
so that ar — 1 AU would always correspond to the distance where the planet receives as much fiux from 
its parent star as the Earth from the Sun. Computations have been done in the range = 0.03— 1 AU. 

Given the star characteristics: mass M^,, radius i?*, effective temperature T^, luminosity L*, spectral 
type Sp and magnitude my, we compute: 

• the orbital radius a in AU; 

• the revolution period P in days; 

• the probability of seeing 3 transits in 150 days: pg x (P < 50 d); 

• the transit duration tr (assuming a mean impact parameter of 0.5); 

• the number of photo-electrons Nph received during tr. 

Then, the SNR detection criterion imposes a minimum value for the signal e = (S/N)/ y/Nph., and 
consequently for the detectable planetary radius Rp — Ri, y/e. 

3.3 Results 

We have plotted the number of detections as a function of the reduced orbital distance for various 
planetary radii (fig.^) and for a false alarm rate of 10~^. For every curve, it is the minimum detectable 
radius (expressed in Earth unit) that is considered. Results are given for the whole mission. This way 
of presenting our results was chosen because of the unknown frequency of the different planetary types. 
Anyone can apply to our curves the scaling factor of his choice. For instance, assuming a 2% probability 
of existence, COROT could detect several tens of "hot Jupiters" enhancing significantly the statistics on 
this class of planets. What is more, COROT has the potential to spot "hot Earths" if any, since for 
example around 4 events would be expected if 20% of the stars exhibit earth planets at 0.05 AU. 
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Raw data 




Dirac comb 



Averaged data 



Figure 1: Raw data are averaged on the duration of the presumed transit, then cross-correlated with a 
Dirac Comb. 




Figure 2: The probabiUty distribution of the cross-correlation products Ck is centered around zero, 
because the absence of a transit is the general rule. A large value of a Ck with respect to ac is indicative 
of a transit with a level of confidence ak- This confidence is also measured by /3, which is supposed 
common to all k. 
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Figure 3: /3 as a function of the false alarm percentage 1 — a. This quasi logarithmic relation allows to 
get high confidence levels with reasonable SNR. 
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Figure 4: This plot shows the number of expected detections during the full life-time of the mission, 
provided 100% of the stars have a planetary companion of radius Rp at the distance a^, and assuming a 
false alarm rate of 10~^. The number attached to each curve is the minimum value of Rp (expressed in 
Earth unit) enabling the detection. 
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